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ABSTRACT 



A set of cosmological variables, which we shall refer to as "supercomoving variables," are 
presented which are an alternative to the standard comoving variables, particularly useful for 
describing the gas dynamics of cosmic structure formation. For ideal gas with a ratio of specific 
heats 7 = 5/3, the supercomoving position, velocity, and thermodynamic properties (i.e. density, 
temperature, and pressure) of matter are constant in time in a uniform, isotropic, adiabatically 
expanding universe. Expressed in terms of these supercomoving variables, the nonrelativistic, 
cosmological fluid conservation equations of the Newtonian approximation and the Poisson equa- 
tion closely resemble their noncosmological counterparts. This makes it possible to generalize 
noncosmological results and techniques to address problems involving departures from uniform, 
adiabatic Hubble expansion in a straightforward way, for a wide range of cosmological models. 
These variables were initially introduced by Shandarin (1980) to describe structure formation 
in matter-dominated models. In this paper, we generalize supercomoving variables to models 
with a uniform contribution to the energy density corresponding to a nonzero cosmological con- 
stant, domain walls, cosmic strings, a nonclumping form of nonrelativistic matter (e.g. massive 
neutrinos in the presence of primordial density fluctuations of small wavelength), or a radia- 
tion background. Each model is characterized by the value of the density parameter f^o of the 
non-relativistic matter component in which density fluctuation is possible, and the density param- 
eter fixo of the additional, nonclumping component. For each type of nonclumping background, 
we identify families within which different values of Qq and flxo lead to fluid equations and 
solutions in supercomoving variables which are independent of the cosmological parameters fio 
and Clxo- We also generalize the description to include the effects of nonadiabatic processes such 
as heating, radiative cooling, thermal conduction and viscosity, as well as magnetic fields in the 
MHD approximation. 

As an illustration, we describe three familiar cosmological problems in supercomoving vari- 
ables: the growth of linear density fluctuations, the nonlinear collapse of a one-dimensional 
plane-wave density fluctuation leading to pancake formation, and the well-known Zel'dovich ap- 
proximation for extrapolating the linear growth of density fluctuations in three dimensions to the 
nonlinear stage. 



Subject headings: cosmology: theory — dark matter — galaxies: intergalactic medium — hydro- 
dynamics — large-scale structure of universe 
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1. INTRODUCTION 



The modern search for an explanation for the origin of the observed structure in the universe has led to 
the need to solve multi-scale, highly nonlinear problems of great complexity which involve the coupling of 
gravitational dynamics, gas dynamics, and cosmological expansion. In what follows, we shall describe a useful 
tool - a change of variables — which simplifies the application of the more familiar results and techniques 
developed to describe a noncosmological gas in a nonexpanding background to problems in cosmology. 

The description of the formation and evolution of large-scale structure in the universe is greatly simplified 
if we restrict our attention to systems with length scales much larger than the Schwarzschild radius of the 
largest mass concentrations in the universe and much smaller than the horizon, c/H, where H is the Hubble 
constant, with nonrelativistic peculiar motions and temperatures. In this limit, known as the Newtonian 
approximation, we can use the ordinary noncosmological, nonrelativistic fluid conservation equations and 
the Poisson equation, modified to include the effects of universal, adiabatic, Hubble expansion, to describe 
the evolution of the baryonic and dark matter that form the galaxies and large-scale structure we observe 
today. 

A straightforward change of variables is sometimes introduced in which the uniform, isotropic Hubble 
expansion is factored out so as to obtain a set of equations that describes the evolution of departures from 
a uniform, isotropic, structureless universe. This is achieved by rewriting the fluid conservation equations 
in comoving variables. These variables are defined by 



and 

Pco-mov = pa{tf , (2) 



where r and p are the proper distance and the mass density, respectively, x is the comoving position, p co -mov 
is the comoving density, and a(t) is the Robertson- Walker scale factor, whose time-evolution is described 
by the Friedmann equation (see eq. [4] below). In these comoving variables, in the absence of structure, 
the mass points are at rest and the density of ordinary, nonrelativistic matter is constant in time. For a 
detailed description of the Newtonian dynamics in comoving variables, see Weinberg (1972) and Peebles 
(1980, hereafter P80). 

The fluid conservation equations in comoving variables contain several terms which distinguish them 
from their noncosmological counterparts. In particular, the momentum and energy equations contain "drag 
terms" that make the peculiar velocities and the thermal energies decay with time, even in the absence of 
departures from uniformity involving pressure or gravitational potential gradients to act as driving forces. 
This reflects the fact that universal expansion, itself, causes adiabatic cooling (i.e. PdV work is done), while 
the kinematics of universal expansion causes a test particle which moves relative to the universal expansion, 
with no peculiar forces acting on it, to appear to fall behind relative to the expanding coordinate system. 

A further change of variables was introduced by Shandarin (1980, hereafter S80), referred to as "tilde 
variables" when modified and applied by Shapiro & collaborators (e.g. Shapiro, Struck-Marcell, & Mclott 
1983; Shapiro & Struck-Marcell 1985; Shapiro ct al. 1996; Valinia et al. 1996). These new variables 
extended the concept of "comoving" to the limit in which the position, velocity, and thermodynamic prop- 
erties (i.e. density, temperature, and pressure) of matter in a uniform, isotropic, adiabatically expanding, 
matter-dominated Friedmann universe all remain stationary. In addition, by replacing proper time by a new 
time variable, the so-called "drag terms" mentioned above disappeared from the momentum equation when 
expressed in terms of these "tilde variable." For an ideal gas with a ratio of specific heats 7 = 5/3, the 
"drag term" in the energy equation disappeared, as well. For such a gas, the complete set of three fluid 
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conservation equations, for matter, momentum, and energy, when expressed in tilde variables, are identi- 
cal to the standard, noncosmological fluid equations. Henceforth, we shall refer to these tilde variables as 
"supercomoving variables." instead, in order to give a more descriptive name. 

The use of supercomoving variables makes it possible to take the known solutions of non-cosmological 
problems and generalize them in a straight-forward way to cosmological situations. An example of this is 
provided by Voit (1996), who used supercomoving variables to address the problem of the evolution of inter- 
galactic blast waves in a matter-dominated universe. The similarity of the fluid equations in supercomoving 
and noncosmological variables also means that analytical methods and numerical simulation algorithms 
which are suitable for studying non-cosmological problems can often be applied to cosmological ones with 
essentially no modification. Examples of the latter include the ID, Lagrangian hydrodynamics method de- 
veloped by Shapiro et al. (1983) and Shapiro & Struck-Marcell (1985) to study cosmological pancakes and 
the new anisotropic version of Smoothed Particle Hydrodynamics (SPH), called Adaptive SPH (ASPH), 
developed recently by Shapiro et al. (1996) to study both cosmological and noncosmological gas dynamics 
by numerical simulation in 2D and 3D. 

Shandarin (S80) also showed that, with an appropriate normalization of the scale factor a(t), it was 
possible to eliminate all dependences, implicit or explicit, upon the density parameter Qq. All matter- 
dominated cosmological models can be grouped into three distinct families, f^o < 1, f^o — 1 > an d > 1, and 
within each family, the fluid conservation equations and their solutions, and the solutions of the Friedmann 
equation, are identical. There is no distinction in supercomoving variables between a model with £Iq = 0.2 
and flo — 0.99, for instance. 

The supercomoving variables introduced by Shandarin (S80) were restricted to matter-dominated models 
(i.e. models in which the mean energy density is dominated by nonrelativistic matter [baryonic and dark], 
so rijco = 0). There are several reasons, however, for considering cosmological models that contain more 
than just the nonrelativistic component. For the particular case of a nonzero cosmological constant, these 
reasons have been summarized by Carroll, Press, & Turner (1992), Ostriker & Steinhardt (1995), and 
Krauss & Turner (1995). The case of a more general, smooth background component has been discussed 
recently by Turner & White (1997). Some of the arguments are as follows. Dynamical estimates of the 
density parameter flo have consistently given values smaller than unity (see, for example, Peebles 1993 and 
references therein) (see, however, Dekel et al. 1993, and references therein for an exception to this trend), in 
contradiction with the most common inflationary scenario (Guth 1981; Linde 1982; Albrctch & Steinhardt 
1982) which requires Oo = 1- (We note, however, that inflationary models with Qq < 1 arc now being 
considered [Lyth & Steward 1990; Ratra & Peebles 1994].) To reconcile the dynamical estimates of fio with 
standard inflation, we must postulate the existence of a second component which must be uniform in space in 
order to escape detection by dynamical methods (hence, we are not referring to the dark matter component 
in galaxy clusters and galactic halos; this component is already included in dynamical estimates of fio)- 
The presence of this uniform component such as domain walls or a cosmological constant also increases the 
age of the universe for a given value of Hq. This helps to reconcile current estimates of Hq, which favor 
Hq « 60-80kms _1 Mpc _1 (Riess, Kirshner, & Press 1995; Freedman 1996; Giovanelli et al 1996), with lower 
limits to the age of the universe based on estimates of globular clusters ages, the well-known age problem. 
With the detection of temperature anisotropy by the COBE satellite DMR experiment, it has become 
possible to fix the amplitude and constrain the wavelength dependence of the primordial density fluctuations 
responsible for large-scale structure (LSS) formation, at least at the superhorizon scale at recombination to 
which the COBE measurement is sensitive (see Bunn & White 1997, and references therein). The much- 
studied, standard, flat, matter-dominated Cold Dark Matter (CDM) model for the origin of galaxies and 
LSS is, however, incompatible with these COBE constraints since it tends to exaggerate the structure on 
the smaller scales measured by the statistical clustering properties of galaxies and LSS in comparison with 
galaxy data (see Ostriker 1993 [§3.2], and references therein; Ostriker & Steinhardt 1995, and references 
therein). The flat CDM model can be reconciled with these COBE results, however, if there is a uniform 
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background component like a cosmological constant which dominates at late times over the nonrelativistic 
matter component. 

The main goal of this paper, therefore, is to generalize the original supercomoving variables of Shan- 
darin, which are limited to matter-dominated models, to include models containing an additional, uniform 
component. In the process, we will demonstrate that different cosmological models corresponding to dif- 
ferent combinations of the fundamental cosmological parameters can be grouped into "families" in which 
the evolution of physical quantities is the same when expressed in supercomoving variables. We will follow 
Shapiro ct al. (1983) and Shapiro & Struck-Marcell (1985) in the inclusion of non-adiabatic processes like 
external heating, radiative cooling, and thermal conduction, and extend their description to include viscosity, 
vorticity, and the effects of magnetic fields in the MHD approximation. 

The remainder of the paper is organized as follows. In §2, we present the fluid equations in noncos- 
mologial variables. In §3, we define the supercomoving variables. In §4, we derive the supercomoving form 
of the fluid equations. In §5, we present the supercomoving form of the collisionlcss Boltzman equation. 
In §6, we derive the cosmic energy equation in supercomoving variables. In §7, we solve the Friedmann 
equation for the various models considered in this paper, and derive the conditions required in order to 
have families of solutions which are independent of the particular values of the cosmological parameters 
for different combinations of these parameters. In §8, we present the solutions of the linear perturbation 
equations for the various models considered. In §9, we discuss exact, nonlinear solutions for the collapse of 
plane-wave density fluctuations, prior to shock and caustic formation. In §10, we present the supercomoving 
form of the Zel'dovich approximation. Summary and conclusions are presented in §11. Our Appendix A 
discusses additional physical processes, including viscosity, heating, cooling, thermal conduction, vorticity, 
and magnetic fields in the MHD approximation. 



2. BASIC EQUATIONS IN NONCOSMOLOGICAL COORDINATES 

2.1. The Cosmological Background 

The evolution of the cosmological background is determined by the relative contribution by its con- 
stituent components to the mean mass-energy density of the universe, as it appears in the Einstein field 
equations. For an isotropic, homogeneous universe with a given composition, this evolution is expressed in 
terms of the Friedmann-Robertson- Walker metric with a cosmic scale factor a(t) whose time-evolution is 
given by the Friedmann equation. 

In this paper, we shall consider two-component models. The first component, called the nonrelativistic 
(NR) component, comprises all forms of matter, luminous or dark, baryonic or non-baryonic, that can 
cluster under the action of gravity, and whose mean energy density varies with time as a(t)~ 3 . The second 
component, which we shall refer to as the X component, does not clump and, instead, has a uniform energy 
density px which varies with expansion according to 

p x (t)<xa(t)- n , (3) 

where n is a non-negative constant. This equation was first introduced by Fry (1985) (see also Charlton & 
Turner 1987; Silveira & Waga 1994; Martel 1995; Dodelson, Gates, & Turner 1996; Turner & White 1997). 
By appropriate choice of n, we recover models with a nonzero cosmological constant (n = 0), domain walls 
(n = 1), string networks (n = 2), vacuum stress (n = 3), or a radiation background (n = 4). In this paper, 
we shall derive a form of the supercomoving variables applicable to all of these various models. Notice that 
models with massive, nonrelativistic neutrinos are also described by equation (3) with n = 3, as long as the 
length scale of the problem studied is much shorter than the damping scale of the neutrinos. For n^3, the 
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X component has a relativistic pressure Px = vpxc 2 , where c is the speed of light and v = n/3 — 1, which 
enters in the Newtonian form of the Poisson equation as an additional source term (see eq. [8] below). 

The time-evolution of the scale factor a{t) is described by the Friedmann equation. For the two- 
component models considered in this paper, this equation takes the form 



= H{tf = H 2 



(1 - ft - Vxo) 



«0 



«() 



+ n 



xo 



«() 



(4) 



(Martel 1995) where we use subscripts to indicate the present value of time-dependent quantities. The 
first term in the right hand side is the curvature term — k/a 2 , where k — — H^a^l — £Iq — flxa)- The other 
two terms represent the nonrelativistic and X component, respectively. The density parameters and fix 
are defined by fl — p/p c and fix = Px / Pc, where p and px = px are the mean density of the nonrelativistic 
and X component, respectively, and p c = 3H 2 /8irG is the critical density, defined as the mean density of a 
flat universe, in the absence of the X component. 



2.2. Fluid Conservation Equations in the Newtonian Approximation 

The dynamics of a self-gravitating fluid in the Newtonian approximation is described by the basic 
hydrodynamical equations for the conservation of mass, momentum, and energy, coupled to the Poisson 
equation and the equation of state. The Newtonian approximation and its validity in cosmology is discussed 
by P80. It is valid in the limit in which all length scales are much smaller than the horizon size, c/H(t), and 
much larger than the Schwarzschild radii of objects within the region of study, and temperatures, equations 
of state, and peculiar velocities are nonrelativistic. In proper coordinates and in the absence of viscosity 
or dissipative processes like thermal conduction or radiative cooling, the continuity, momentum, and energy 
equations, and the equation of state take the form 

|? + V-(pu)=0, (5) 
^ + (u-V)u = -V$-^, (6) 

f| +u . Ve = -^V.u, (7) 
at p 

p=(j-l)ep, (8) 



where p is the matter density, e is the specific internal energy, p is the gas pressure, u is the velocity, $ is 
the gravitational potential, and 7 is the ratio of specific heats. These equations also assume that there is no 
interaction besides gravity between the nonrelativistic component and the X component. In Appendix A, 
we generalize these equations to include additional physical processes, such as heating, cooling, thermal 
conduction, viscosity, and magnetic fields. 



2.3. Gravitational Potential in the Newtonian Approximation: Modified Poisson Equation 

The Newtonian form of the Poisson equation must be modified to account for the presence of the 
X component. To derive the Poisson equation, we consider the acceleration g a between two freely- falling 
particles separated by a proper distance £ Q . In the small-scale limit, where £ a is much smaller than the radius 
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of curvature of the universe, the space-time is locally Minkowskian, and the equation for the acceleration 
reduces to 

9 a = R a oo/sf , (9) 

(Peebles 1993, p. 268) where R l jki is the Riemann-Christoffel curvature tensor, and where we sum over 
repeated indicies. By taking the divergence of equation (9), we get 



V 2 $ 



-V • e = AttGT° 



(10) 



where T is the stress-energy tensor. For the models we consider, the total stress-energy tensor in the fluid 
rest frame is given by 

-p + Px " 

px/c 2 
px/c 2 

. Px/c 2 . 



(11) 



where px and px are the energy density and pressure of the X component, respectively, and c is the speed 
of light. In equation (11), we neglected the nonrelativistic pressure term p/c 2 , which is smaller than p by a 
factor of order (c s /c) 2 , where c s is the sound speed. The evolution of the density px of the X component is 
given by 

Ta (Pxa ) = -— (12) 
(Weinberg 1972, p. 472). By combining equations (3) and (12), then yield 



Px n 



(13) 



Equations (10), (11), and (13) reduces to 

V 2 $ = 4ttG p+(n-2)p x 



(14) 



This is the final form of the Poisson equation in the presence of the X component. Notice that the contribution 
of the X component to the gravitational source term is negative for n < 2. This implies that the gravitational 
force resulting from the presence of domain walls or a cosmological constant is repulsive. 



3. SUPERCOMOVING VARIABLES 

One of the most desirable features of the comoving and supercomoving variables is that in these variables, 
in the absence of structure, the density is constant in time and mass elements are at rest. In supercomoving 
variables, a gas with ratio of specific heats 7 = 5/3 also has its thermodynamic variables remain constant in 
the absence of structure. We preserve these two properties when we generalize supercomoving variables to 
the cases involving an additional nonclumping background component, by making the following definitions, 



a 3 p 
p* 



(15) 
(16) 
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where the peculiar velocity v is defined by 



dt = 



av 

V* 

dt 
a 2 t* 
a 2 ^ 



5 

a p 



a 2 e 

£* 



u-JJr. 



where H = a 1 da /dt is the Hubble constant, and the peculiar gravitational potential <fi is related 
Eulcrian gravitational potential $ by 

^_ 2nGpr 2 | 2(n-2)^r 2 | SlgH* / ag\ 3 ^ | »^ ff 2 (n - 2) /ap^ | 

3 3 4 \ a / 4 \ a 



(17) 
(18) 
(19) 
(20) 
(21) 

(22) 
to the 

(23) 



This is a standard transformation, except for the second term which we introduce to take into account the 
presence of the X component. In equations (15)— (21), the quantities r„, p*, v„, t*, cj>*, p*, and are fiducial 
values. Only is independent. The other fiducial quantities are defined by 



p* = 

— 

v* = 



P* 



8ttG ' 
2 / /„ ^ 



r* 
t * 

2 / ,2 2 

P*rl/tl = p*v 2 , 
= vl , 



(24) 

(25) 

(26) 

(27) 
(28) 
(29) 



where po is the cosmic mean matter density at the present epoch and the parameter /„ is defined by 

1, n^3; 

n = 3. 



fn 



fin 



(30) 



fid + fixo 



For any particular cosmological model, we can eliminate a(t) using the Friedmann equation, and inte- 
grate equation (18) to get an explicit relation between t and t. For future reference, equation (18), together 
with the Friedman equation (4), implies the following expressions for the Hubble constant and the first and 
second derivatives of the scale factor, 



H(t) 

'da 
. dt 



1 da - 1 H 

a 3 t* dt a 2 t st 



.2 tt2K 

t * tL o a o 



(1 - fio - Vxo) 



fir. 



fi 



XQ 



«0 



6— n 



(31) 
(32) 
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2(l-fi -fixo) 



351o / a 
~2~ W 



2/ V»o 



5— n 



(33) 



where equation (31) defines the quantity H = a 1 da/dt, the supercomoving Hubble constant By combining 
equations (31) and (32), we can eliminate da/dt and express H as a function of a. 



4. SUPERCOMOVING FLUID EQUATIONS 



4.1. Transformation of Derivatives 



To convert the hydrodynamical equations (5)-(8) and (14) to supercomoving variables, we need to 
express the derivatives relative to r at fixed t and t at fixed r as functions of the derivatives relative to f at 
fixed t and t at fixed f , using equations (15) and (18). This yields 



df 
Of 



(V/) t = 



1 



1 



df 
dt 



(v/) t -, 



Wr-(V/) t - 



(34) 
(35) 



where V is the gradient relative to f . 



4.2. The Continuity Equation 

To derive the form of the continuity equation in supercomoving variables, we substitute equations (15), 
(16), (17), (22), (24), (26), (34), and (35) into equation (5). After some algebra, we get 



(36) 



This equation has exactly the same form as equation (5). 



4.3. The Momentum Equation 

We substitute equations (15), (16), (17), (19), (20), (22)-(28), (34), and (35) into equation (6), and 
eliminate the terms in da/dt and d 2 a/dt 2 using equations (32) and (33). We get, after some algebra, 



f + (v.V)v = -?-W. 
at p 



(37) 



Again, we have obtained an equation that has the same form in supercomoving variables as in noncosmological 
variables. The corresponding equation in comoving variables contains an additional drag term. 
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4.4. The Poisson Equation 

We substitute equations (16), (19), (23), (24), (25), (27), and (35) into equation (14), and get 

V 2 ^6a/ rl (|-l) , (38) 

where we used p = a?pj p* = a 3 p/p = Oq (sec cqs. [16] and [23]) to eliminate a in equation (24). For n = 3, 
we can rewrite equation (38) as 

V 2 = 6a(^-lV (39) 

VPtot / 

where p to t — P + Px is the total density, including the X component. The reason for treating the case n = 3 
differently will be explained in §7.3.3 below. 



4.5. The Equation of State 

We substitute equations (16), (20), (21), (24), (28), and (29) into equation (8). The equation of state 
becomes 

p=(7-l)/5e, (40) 

which has the same form as equation (8). 



4.6. The Energy Equation 

We substitute equations (15), (16), (17), (20), (21), (22), (24), (26), (28), (29), (34), and (35) into 
equation (7). The energy equation becomes 

-=+v- Ve + W(3 7 - 5)e = -^V-v. (41) 
at p 

This equation contains a drag term, which vanishes only for 7 = 5/3. We can derive the evolution of the mean 
specific internal energy e by considering a universe with no fluctuation. In this case, v = and d/dt = d/dt. 
Equation (41) then becomes 

d~e 

-=+W(37- 5)1=0. (42) 
dt 

The solution is 

\ 5-37 
a N 



m = ~^-) . (43) 

Hence, the mean specific internal energy decreases with expansion for 7 > 5/3, and increases for 7 < 5/3. 
Since the mean mass density does not vary with expansion in supcrcomoving variables, the equation of state 
gives 

/ \ 5-37 



#*)=*(-) (44) 
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The reason for the absence of drag term in equation (41) for the case 7 = 5/3 only is easily understood. 
The momentum equation (37) contains no drag term. From this equation, we can derive an equation for the 
evolution of the specific kinetic energy v 2 /2, which contains no drag term either. In the case of 7 = 5/3, the 
internal energy of the gas is nothing more than the microscopic specific kinetic energy of the particles that 
compose the gas. Hence, because of the absence of drag term, the mean specific internal energy remains 
constant as the universe expands. But if 7 is not 5/3, then there is more to the internal energy of the gas 
than just the specific kinetic energy of the particles. For instance, for diatomic gases with 7 = 7/5, some of 
the internal energy is in the form of molecular rotation. Hence, the drag term in the energy equation (41) 
represents the component of the internal energy that is not associated with the microscopic motion of the 
atoms and molecules in the gas. 



5. THE SUPERCOMOVING COLLISIONLESS BOLTZMANN EQUATION 

The equations derived in §4 describe the evolution of a collisional gas. These equations cannot be used 
for describing the evolution of a collisionless system of particles, because the velocity u is a multivalued 
quantity. The proper way to describe such system is in terms of a distribution function /(r, u, t) which 
is defined as the number density of particles in phase space. The evolution of the distribution function is 
described by the collisionless Boltzmann equation, also known as the Vlasov equation. In noncosmological 
variables, this equation takes the form 

df 

^+u-V/-V$-V u / = 0, (45) 

where the gravitational potential is obtained by solving the Poisson equation, which, for collisionless systems, 
is given by 

V 2 $ = 4irG 



J fd 3 u 



+ Pb 



(46) 



where we included the baryon density pb to take into account systems in which the nonrelativistic component 
is a mixture of a collisionless and a collisional (gaseous) component. After switching to supercomoving 
variables, the Boltzmann equation becomes 

^+v-V/-W-Vv/ = 0, (47) 
where /(f , v, i ) = /(r, u, t), and the Poisson equation becomes 




P**+h -lV n^3: 



JfcPv + px+Pb 
(Jfd 3 v)+p x + ~p b 



(48) 



- 1 , n = 3; 



where ( ) indicates space-averaging. 



6. THE COSMIC ENERGY EQUATION 



Using the fluid equations derived in the previous sections, we can derive an equation for the evolution 
of the total energy. Our derivation follows the same lines as the one presented in P80 for the cosmic energy 
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equation in comoving variables. Taking the dot product of pv with the momentum equation (37), we get 



10(^) ^P + v. (/5v) 



2 dt 



at 



V • [pv + p<t>v + 



pv 2 v 



- pV • v + </>V • (pv) . 



(49) 



Using the continuity equation (36), we cancel the term in square brackets, and eliminate V • (pv) in the 
last term. We also use the energy equation (41) to eliminate the pV • v term. After some additional vector 
manipulation, we get 



1 d(pv 2 ) 

2 dt 



V • [ pv + p(\><f 



pv 2 v 



pev 



djpi) 
dt 



+ e 



- H(3 7 - 5) pi - ~4>% . (50) 
ot 



Again, we cancel the term in square brackets using the continuity equation. We then integrate equation (50) 
over all space. For a system with finite extent, the divergence term vanishes. We now introduce the following 
definitions for the total kinetic energy and total thermal energy of the system: 



E = 



ped 3 r . 



(51) 
(52) 



With these definitions, equation (50) reduces to 

dt dt 'J 



dt 



(53) 



Notice that integration over all space has transformed the partial time derivatives into total time deriva- 
tives, and that we introduced the constant mean density p in the last term. To compute the last term in 
equation (53), we first solve the Poisson equation (38), and get 



4>(v 



, = 3afn f 
' 2tt~p J 



r — r 



(54) 



We can easily show that this is indeed the correct solution by directly comparing equations (38) and (54) to 
their Eulerian counterparts. With this solution, the last term in equation (53) becomes 



d (f>-~p) d 3~_ ad 
dt 2 dt 



fi{f)[p{f)-p\<Pi 



(55) 



We then introduce the following definition for the gravitational potential energy of the system: 

U=l[ 4>(P~~P)d 3 r. 



(56) 



Equation (53) reduces to 



—,{f + E + U)=H (5 - 3 7 )£ + U 
dt 



(57) 



This is the final form of the cosmic energy equation in supercomoving variables. The corresponding equation 
in comoving variables (see P80, eq. [24.19]) is significantly more complicated. 
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7. SOLUTIONS OF THE FRIEDMANN EQUATION 



7.1. Cases with Elementary Solutions 

The evolution of the cosmological background is described by the Friedmann equation (4), written in 
supercomoving variables in equation (32). For convenience, we rewrite equation (32) as follows, 

2 = 4(1 -"T^ + 4,3 + l^V- , (58) 
where 

x = % a , (59) 

t = (a f n ) 1/2 t. (60) 

The solution of this equation is 

1 f , /1-Oo-Oxo 4 . 3 . Qxo 6-n\ ^ i a -\\ 

T= 2j dx { n x+x+ ^ x ) ■ (61) 

This integral has no elementary solution for arbitrary values of the parameters Oo, ^xo, an d n. There are, 
however, several physically interesting cases for which elementary solutions exists. We can identify the cases 
that have or might have an elementary solution by studying the properties of equation (61). 

First, all known cases which are physically relevant have < n < 4, with n integer. Hence, we can 
rewrite equation (61) as 

r=U -P{x)-^\ (62) 



2 J x 
where 

P{x) = i-no-n^ + x + nxo x ,- n (63) 

is a polynomial of degree 4 or less. If n = 2, 3, or 4, P(x) is a second degree polynomial, and equation (62) 
always has an elementary solution. If n = 1 (domain walls) or n — (cosmological constant), then P(x) 
is a polynomial of degree 3 or 4, the integral in equation (62) is elliptic, and does not have an elementary 
solution in general. There are special cases, however. Suppose that P(x) has a double root at some particular 
value x — x st . We can then rewrite P(x) as 

P(x) = (x - x st ) 2 Q(x) , (64) 

where Q(x) is a linear polynomial if n = 1 or a quadratic polynomial if n = 0. Substituting equation (64) 
into equation (62), we get 



X\X - X st I 



for which there is always an elementary solution. The necessary condition for the polynomial P{x) to have 
a double root is that P(x) and its first derivative both vanish at which becomes 



1 — — ^Jfo\ 2 , , ^xo 4 

n K + ' Tst + -ri ' 



x 2 st+ x st + -^x 4 st - n = 0, (66) 



^ st + l + (4-n)-^.x s 3 t -" = 0. (67) 
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We can eliminate x 8t from these equations to get a condition relating fl and £lxo- We get, after some 
algebra, 



1 — f2o — Qxo 
The value of x s t is then given by 

Xgt = < 



/^ao^oV^ 3 n 
-3| : — - I , n = 0; 



(68) 



-2(n X0 n ) 1/2 , n=l. 



go 
2ft 



A"0 



1/3 



1/2 



, n = 0; 



n = 1 



(69) 



The physical meaning of x s t is easily understood. By combining equations (58), (66), and (67), we can show 
that a and a both vanish at x = x st ; hence the universe is static at a = a x st . The case n = is known 
as the Einstein static model. In this model, the repulsive force caused by the cosmological constant exactly 
cancels the gravitational force. As we see, the possibility of a repulsive force canceling gravity also exists 
for domain walls (this is also shown by eq. [14], in which the X component acts as a negative density for 
n < 2 models). Since the universe is clearly expanding at present, x must be smaller than the equilibrium 
value x s t. Equations (58) and (68) then describe a universe which decelerates and asymptotically approaches 
a static universe as x gets closer to x st . 



7.2. Families of Solutions 

Equation (61) (or [32]) relates the cosmic scale factor a to our supercomoving time variable t in any 
particular model. According to this equation, models with different values of fl and Qxo can often be 
grouped into "families" within which the dependence of a on i is the same, as follows. 

There is a basic difference which occurs between those cases for which P(x) defined by equation (63) 
reduces to only two terms and those cases for which P(x) has three terms. Let us call the former case, 
category 1, while the latter is category 2. For cases in category 1, it is always possible to find solutions 
with no explicit dependences upon the cosmological parameters fto and Qxo- If P( x ) contains only two 
terms, the change of variable x — by, with an appropriate choice for b, will make the coefficients of these 
two terms identical. We can then move that coefficient to the left hand side of equation (62) and eliminate 
it by rescaling r. Cases for which P(x) has only two terms include matter-dominated models (flxo — ; as 
described in S80), any zero-curvature model (1 — fto — ^xo = 0) regardless of the presence of an additional 
uniform background component, models with cosmic strings (n = 2), and models with a nonrelativistic 
uniform component (n = 3). In each of these cases in category 1, the solution for a(t) as a function of t 
reduces to a finite number of functions which describe the infinite number of combinations of values for VLq 
and Slxo which can occur in each case. We present these solutions in §7.3. In all other cases, P(x) contains 
3 terms (category 2), and there are an infinite number of solutions for the dependence of a(t ) on t, but we 
can still identify combinations of f^o and &>xo which share a common solution. To find these combinations, 
we set x = (flo/flxo) 1 ^ 3 ^ n ^y m equation (63), which makes the coefficients of the last two terms identical. 
We get 

>(n-2)/(3-n) 



P(y) 



«o \ 1/(3 - n) \ (i - n - nxoM 1 -"'^-"' 2 



Q 



AO 



l/(3-n) 
AO 



y~ + y + y 4 ~ 



(70) 
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As long as the following condition is met: 



k = 



(i-n -n X0 )n("- 2)/(3 - n) 



= const 



(71) 



,l/(3-n) 
'XO 



the solution will not depend upon Qq and Qxo, For any combination of f^o an d ^xo which satisfies equa- 
tion (71), there is a single dependence of a(t) on t, which is different for different values of n. There are, 
therefore, an infinite number of such cases, corresponding to the infinite range of values of k. The particular 
cases with n = or 1 for which closed- form, analytical solutions exist (cq. [68]) do satisfy the condition in 
equation (71). 

Henceforth, we shall refer to the combination of values of f^o, Qxo, and n which share a common 
solution for the dependence of a(t ) on t as "families." The matter-dominated cases (Clxo = 0) considered by 
Shandarin (S80) belong to category 1 of the cases described above. As shown by Shandarin, these models 
which are parametrized by £Iq are, in supercomoving variables, fully described by three families, closed, flat, 
and open, and within each family, the dependence of a(i) on i is identical, independent of the value of CIq. 
We have shown above that this concept of "families" can be extended to include the much wider range of 
cosmological models considered here, parametrized by D,q, Qxo, & n d n. For category 1 models, of which the 
matter-dominated are one example, there are also just a finite number of families per model (i.e. for each 
value of n) as in the matter-dominated case. For category 2 models, however, each value of n corresponds 
to a different family. This might look like a mathematical trick at first sight, where we use equation (71) 
to group different models that appear to be unrelated. However, the fact that, within a given family, the 
models will be identical in supercomoving variables indicates that the various models assembled into the 
same family are indeed deeply related. For instance, for any particular value of n, the flat (zero-curvature) 
models constitute a family, defined by n = 0. Hence, all flat models with the same value of n will have 
the same solution for a(t ) in supercomoving variables. Other physically interesting families include critical 
models with a cosmological constant, defined by n — 0, n — — 3/4 1 / 3 , and critical models with domain walls, 
defined by n — 1, n = —2. Now that we have derived equation (71) and defined the concept of family, we 
are ready to solve the Friedmann equation for the various cases which have elementary solutions. 



7.3. 



Solutions 



7.3.1. Matter- Dominated Universe (ilxo = 0) 



These are the only models considered in S80. The solutions of the Friedmann equation are 




(t 2 - 1) 
t- 2 , 
(P + l) 



-l 



l 



fio < l; 
Oo = l; 
fi > 1, 



(72) 



where the present value a of the scale factor is 




(73) 



These values of oo were chosen in order to eliminate the dependences upon f2 in equation (72). These 
supercomoving- variable solutions are remarkably simple compared with their comoving- variable counterparts. 
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In equation (72), i is negative and increases with time from the value i = — oo at the Big Bang. For the 
fiat (Oo = 1) model, t = — 1 at the present, and approaches as the universe expands to infinite radius. 
For the open (O < 1) model, t < — 1 at the present, and approaches —1 as the universe expands to infinite 
radius. For the closed (Oo > 1) model, t < at present, the universe reaches a maximum expansion at t = 
and recollapses at t = oo ("the Big Crunch"). The solutions of equation (72) for the cosmic scale factor 
a(t) as a function of t are plotted in Figure la (solid curves). The solid dot on the curve for the flat model 
indicates the present. For all other models, the location of the present is a function of Oo. In Figure lb, we 
have plotted the "supermoving Hubble parameter" H = a~ 1 da/dt defined by equation (31), also using solid 
curves. 



7.3.2. Universe with Infinite Strings (n = 2) 

This case is extremely simple. Setting n = 2 in equation (61) has exactly the same effect as setting 
Oxo = 0. Hence, the solutions for a universe with strings are exactly the same as equations (72) and (73) 
above for a matter-dominated universe, independent of the value of Oxo- This results from the fact that a 
uniform background term with n = 2 behaves exactly like the curvature term in the Friedmann equation. 
Hence, the solid curves plotted in Figure 1 for the matter-dominated models apply to these models as well. 



7.3.3. Universe with a Nonclumping, Nonrelativistic Component (n = 3) 
This case is also very simple. We define an effective density parameter O as 

O = O + Oxo - (74) 



Equation (61) then reduces to 



"(^^ ■ <75) 



This equation has exactly the same form as in the case of a matter-dominated universe, except for the extra 
factor (Oo/Oq) 1 / 2 in front of the integral. However, this extra factor disappears after we eliminate r using 
equations (30) and (60), and we get exactly the same solutions as for the matter dominated case, 




O + Oxo < 1 
Oo + Oxo = 1 
O + Oxo > 1 



(76) 



where the present value of the scale factor is given by 

1 — Oo — Oxo 



0„ + Oxo ' n ° + n *°< 1 ; 
a = < 1 , O + Oxo = 1 ; (77) 



Oo + Oxo- 1 
Oo + Oxo 



O + Oxo > 1 • 



This was the reason for treating the case n = 3 differently in equations (25) and (60). Otherwise, the 
solutions would have shown dependences upon O and Oxo- The solutions (76) and (77) are identical to 
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the solutions (72) and (73) for the matter-dominated cases (and hence, the solid curves plotted in Figure 1 
apply to these models as well) if we replace O in the latter by £!(,. This makes sense since, for n = 3, 
the clumping, nonrelativistic and nonclumping, nonrclativistic X components are indistinguishable in the 
Friedmann equation. They are distinguishable in the Poisson equation, however, since one component can 
contribute to the peculiar gravitational potential while the other cannot. Hence, we end up with a different 
form for the Poisson equation, which is the price we must pay for having no dependence in equation (76) 
upon the cosmological parameters. 



This modification of the Poisson equation is actually physically interesting. We can rewrite equation (39) 



as 



|V 2 ^| 



6 a 


P+~Px 


- 1 


< 6a 






P + Px 






p 



(78) 



Since the solution a(t) is the same for both matter dominated and n = 3 models, the interpretation of the 
inequality in equation (78) is that the presence of the X component effectively weakens the peculiar gravity. 
This explains, for example, why the growth of density fluctuations in the baryons and CDM is reduced in a 
universe which includes massive neutrinos, the Cold + Hot Dark Matter (CHDM) model, compared with the 
growth in a CDM model with the same f^o, for wavelengths smaller than the neutrino free-streaming length. 
In the CHDM model, the growth of small-wavelength fluctuations occurs as if the neutrinos corresponded 
to a nonclumping, nonrelativistic background with Q, v = flxo, until very late times when long wavelength 
neutrino density fluctuations achieve significant growth too. 



7.3.4- Universe with a Radiation Background (n = 4) 
The solution of equation (61) for n = 4 is 



l) 2 fijc (l - fi - Vxo) 



-I -1 



fig 



(79) 



where the present value a® of the scale factor is 



a 



go 



xo 



(80) 



As in other cases, this solution gives a = (Big Bang) for t = — oo. The behavior of the solution depends upon 
the curvature of the universe. If flo + flxo < 1, a becomes infinite when the square bracket in equation (79) 
is zero. This gives a maximum value for t, 



t(a 



oo 




4ft X0 (i -no- n xo ) 



-,1/2- 



(81) 



If Clo + Clxo > 1, the square bracket in equation (79) can never be 0, implying that a can never be 
infinite. This case corresponds to that of a bound universe. The Big Crunch (a = 0) occurs at t = +oo. The 
time at turnaround can be found by setting da/dt — 0. We find 



^turnaround 



-im 

4 



1 + 



4fi xo (^o + ttxo - 1) 



(82) 
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Unlike in the previous cases, we have not succeeded in obtaining a solution that has no dependences 
upon the parameters flo and flxo- However, the solution (79) will be the same for different models with the 
same value of k = fixo(l — — ^xo)/^o, which is precisely what equation (71) predicted. The essential 
difference between this case and the previous cases is that matter-dominated models, models with cosmic 
strings, and models with a uniform nonrelativistic component are one-parameter models (the parameter 
being fl for the first two cases, and Sl = Sl + Qxo for the latter case). It is then possible, as we saw, 
to find solutions with no explicit dependence upon that parameter. The radiation model, however, is a 
two-parameter model, and we cannot eliminate these parameters from the solution. The solution to this 
problem, as we saw in §7.2, is to consider models such as the radiation model, not as two-parameter models, 
but as an ensemble of one-parameter models, each model being characterized by the value of the constant k 
in equation (71). Equation (79) then reduces to 

V 2t ~-i) 2 



a = e~ 2t 



— K 



(83) 



For any particular value of k, we are in a family for which the cosmological parameters f2o and Qxo arc 
related by equation (71). All these models have the same solution, equation (83), independent of the value 
of the cosmological parameters. We have plotted the solution for flat models (k = 0) in Figure 1 (dotted 
curves). The radiation component dominates at early times, which explains why this model has a different 
behavior at early time compared to all other models. The most interesting family is the one for which the 
spatial curvature is zero, which corresponds to k = 0. In the limit t — > 0~, we recover the solution a = t~ 2 , 
for the flat matter-dominated model, since the relative contribution of the radiation component decays as 
the universe expands. 



7.3.5. Critical Universe with a Nonzero Cosmological Constant (n = 0, Cixo = ^xo.crit) 

A universe with a positive cosmological constant is marginally bound if f2o and flxo satisfy equation (68) 
for n = 0. We substitute this expression into equation (61). The resulting integral has an elementary solution, 
which is 

/ 1\ 1/2 1 o+l- [3a(a + l)l 1/2 , 
t= 1 + 7 -^77i ln r. 777?' (84) 



a 



3 1/2 a+l+ [3a(« + l)] 1/2 ' 



where the present scale factor is given by 



O \ 1 / 3 



"^l^ 1 ■ < 85) 



Unfortunately, we cannot invert equation (84) to express a as a function of t (except numerically). In this 
model, t = — oo at the Big Bang, and goes to t = +oo as the scale factor a approaches its maximum value 
flmax = 1/2- We have plotted this solution in Figure 1. 



7.3.6. Critical Universe with Domain Walls (n = 1, flxo = Slxo.crit) 

A universe with domain walls is marginally bound if f^o and Qxo satisfy equation (68) for n = 1. We 
substitute this expression in equation (61). The resulting integral has an elementary solution, which is 

t = l-r + tanlT 1 a 1/2 , (86) 
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where the present scale factor is given by 



o \ V2 



Again, we cannot invert equation (86) to express a as a function of t. In this model, i — —00 at the Big 
Bang, and goes to t = +00 as the scale factor a approaches its maximum value a max = 1. We have plotted 
this solution in Figure 1 



7.5.7. Flat Models with Domain Walls or a Cosmological Constant (fio + &xo = 1, n = 1 or 0) 

These cases are interesting and important for several reasons. First, the flatness of the cosmological 
model is usually regarded as a requirement of the standard inflationary scenario. Second, as summarized by 
Ostriker & Steinhardt (1995), the flat model with nonzero cosmological constant and Cold Dark Matter is 
of great interest as an explanation for several observations which are otherwise discordant with the standard 
flat CDM model without a cosmological constant. For example, the presence of a cosmological constant (or 
domain walls) increases the age of the universe for a given observed H , helping to solve the so-called age 
problem. As a practical matter, therefore, it is convenient that, according to the discussion presented in 
§7.2, solutions of the Friedmann equation for flat models, elementary or not, will have no explicit dependence 
upon the cosmological parameters. 

Equation (61) for flat models reduces to 

HM^^T' 2 - (88) 



Using equations (59) and (60), we can write this integral as 

dy 



2 J 1 t/ 3 / 2 (l + y3-n)l/2 ' 



(89) 



where the present scale factor is defined by 



O \ V(3-«) 



«o=(-jf) (90) 



By writing equation (89) as a definite integral, we have fixed the value of the integration constant by 
imposing (arbitrarily) a = 1 at t — 0. As expected, equation (89) has no explicit dependence upon fig 
and flxo- Unfortunately, this integral has no elementary solutions for the cases n — and n = 1, and thus 
has to be evaluated numerically. In the limit a <C 1, the term y 3 ~" can be neglected, and we recover the 
solution t — —a 1 / 2 , a = t~ 2 of the matter-dominated flat universe (eq. [72]). This was expected, since the 
contribution of the uniform component is negligible at early times for n < 3. The value of t at present is 
obtained by integrating equation (89) numerically for a = ao- The results are given in Table 1 for various 
cases. Finally, t approaches a finite value as the universe expands to infinite radius. This value can be 
computed numerically by integrating equation (89) with a = 00. The results are also given in Table 1. 
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7.3.8. Solution Families 



The results presented in this section are summarized in Figure 2, where we plot the relationships between 
fio and Qxo which define family membership for each value of n, for models with n = 0, 1, 2, 3, and 4. 
For models with n = 0, 1, and 4, families are defined by equation (71), with each possible value of k 
corresponding to a particular family. Families are thus represented by curves of constant k. We have plotted 
several of these curves in Figure 2a, 2b, 2c, and 2f, and indicated the corresponding values of k. Thick 
dashed curves indicate particular families that were discussed in this section (fiat and critical models with 
a cosmological constant or domain walls, and flat models with radiation). The cases n = 2 and n = 3 are 
different. As we saw, there are only 3 families for these models. These families are represented in Figure 2d 
and 2e by the dashed line, and the half-planes on cither one side of it. The dashed areas in Figure 2a and 2c 
represent regions in parameter space that are excluded by the existence of the Big Bang. In these regions, 
the Friedmann equation (4) predicts a minimum value a — a m ; n > for the scale factor. For a < a m j n , 
equation (4) gives (a) 2 < 0. 

The solutions presented in this section describe the evolution of the cosmological background. The next 
step consists of describing the growth of fluctuations leading to the formation of large-scale structure. This 
problem is much too complex to be solved analytically for the general case. In the next two sections, we 
consider two particular cases of great cosmological importance for which analytical solutions exists, either 
because various terms in the fluid equations can be neglected, or because the problem has a high degree of 
symmetry. These two problems are: (1) the growth of fluctuations in the linear regime, and (2) the growth 
of a ID, plane-wave density fluctuation from the linear to the nonlinear regime, leading to the formation of 
a collisionless pancake. The extension of these results to the nonlinear regime for a general distribution of 
initial density fluctuations is then made possible by the well-known Zel'dovich approximation presented in 
supercomoving variables is §10. 



We can use the fluid conservations equations in supercomoving variables to describe the growth of small- 
amplitude density fluctuations in the early universe. First, we decompose the density p and the pressure p 
into a uniform and a space- varying component, 



where p and p are the mean values of the density and pressure. In supercomoving variables, p is constant 
and p is a function of time only (see eq. [44]) (and p = constant if 7 = 5/3). We substitute equations (91) 
and (92) into the fluid equations. These equations reduce to 



8. LINEAR PERTURBATION THEORY 



P = ~p(l + S), 
p = p + Sp, 



(91) 
(92) 



+ V.[v(l + *)]=0, 




(93) 



-V0, 



(94) 



V 2 4> = 6af n 6 . 



(95) 



A trivial case which satisfies these equations isc) = O,v = O,0 = O, Sp = 0, which corresponds to a universe 
with no fluctuation. Assuming that the fluctuations are small, we can solve equations (93)-(95), in general, 
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using perturbation theory. After dropping second order terms, equations (93) and (94) reduce to 



OS 
dt ' 
9v 
~dt 



V-v = 0, 
V5p 



- V(f>. 



(96) 
(97) 



(Eq. [95] is already linear.) To solve these equations, we take the divergence of equation (97), then use 
equation (96) to eliminate v, and equation (95) to eliminate (j). We get 



(PS 

dt 2 



-~- = — h 6af n d . 



(98) 



To solve this equation, we need a relationship between the pressure fluctuation dp and density contrast 6. 
For adiabatic perturbations, V 2 <5p = (? s pV 2 8 to first order, where c s = (^p/p)^ 1 ^ 2 is the mean sound speed. 
Equation (98) reduces to 



^ = ty 2 S + 6af n S. 



(99) 



8.1. Critical Jeans Length 

The density contrast 5 can be expressed as a sum of plane waves, 



ik-r 



(100) 



We substitute equation (100) in equation (99). The terms for each value of k must cancel out separately. 
We then get a separate equation for each mode k, 



dt 2 



= 6a/„ 



<5k , 



(101) 



where A = 27r/|k| is the wavelength of the mode, and Xj is the Jeans wavelength, defined by 

Xj = 27r(6a/„)- 1/2 5 s . 



(102) 



Perturbations can grow only if A > Xj. Modes with A < Aj result in sound waves. The mean sound speed is 
constant for 7 = 5/3 models, and a function of a alone for 7 ^ 5/3 models. Thus, the only time dependence 
in equation (102) is contained in the scale factor a(t). This implies that within any given family, the value 
of the Jeans length is independent of the cosmological parameters. For instance, we can compute the Jeans 
length for matter-dominated models by substituting equation (72) into equation (102). We get 



Xj = 



2tt - 



(3 7 -5)/2 



( f ~2 _ 1)37 /2-2 5 Q Q<1 . 

■ < t 3 ^- 4 , fi = 1 ; 

(f -2 + 1)37/2 -2 5 no>l 



(103) 



We note that, if the temperature of the baryon-electron gas remains coupled to that of the cosmic 
microwave radiation background by Compton scattering, then 7 = 4/3 is the appropriate value to use in 
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equation (43) and c s oc a 1 / 2 . In that case, equation (102) indicates that Aj is independent of time. This 
means that the baryon Jeans mass Mj (in physical units, where Mj oc pAj) must be constant in time, as 
well, in that case. This latter result is a well-known one for the Einstein-de Sitter model, and is relevant to 
the postrecombination epoch between z <~ 10 3 and z ~ 10 2 for mean baryon densities consistent with Big 
Bang nucleosynthesis constraints for the observed values of the Hubble constant (see, for example, Shapiro, 
Giroux, & Babul 1994). However, by solving the problem in supercomoving variables, we have generalized 
this result to a very much wider class of models. 



8.2. Zero Pressure Solutions 

In the limit A » Aj, equation (99) reduces to 

8 2 5 



di 2 



Qaf n 5 . 



(104) 



This equation has a slightly simpler form in supercomoving variables than in comoving variables. However, 
both equations are equally complicated to solve. Indeed, the first step in solving equation (104) is to replace 
the independent variable t by x — a/ao using equations (25), (32), and (33). Equation (104) reduces to 



[(i - o - n xa )x 3 + n x 2 + n X0 x 5 - n ]6" 
+ 



2(1 - n - n XQ )x 2 + ^n x + (s - 1 ) n X0 x 4 -' 



5 = — S -> 



(105) 



where a prime stands for d/dx (notice that /„ has disappeared from the equation). This equation is the 
same in supercomoving and comoving variables. It has two independent solutions, a decaying solution 5- 
and a growing solution 5 + . We shall focus on the growing solution, which is responsible for the formation 
of structures. Numerous solutions of equation (105) have been published. Hence, we only need to reexpress 
these known solutions in terms of supercomoving variables, by eliminating the scale factor a(t ) using the 
solutions derived in §7.3. 



8.2.1. Matter-Dominated Universe (Vt xo = 0) 

The solutions for matter-dominated models are given in Weinberg (1972) and Groth & Peebles (1975). 
In terms of supercomoving variables these solutions are 



5+ 



l + 3(t 2 -l) 1-Uanh 



1 + 3(t 2 + 1) 1 ~Uan 



1) 


, ^0 


< 1 










= 1 


i) 


j ^0 


> 1 







(106) 



where the factors of 5/2 were chosen such that all solutions reduce to 5 + = tr 2 at early time. These solutions 
were given in S80 (Actually, S80 eq. [15] contains some typographical errors. We have corrected this in our 
eq. [106] above). 
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8.2.2. Universe with Infinite Strings (n = 2) 

As in the case of the Friedmann equation, setting n — 2 in equation (105) has the same effect as 
setting Clxo = 0, hence the solutions (106) also apply to models with infinite strings. 



8.2.3. Universe with a Nonzero Cosmological Constant (n = 0) 

The general solutions of the perturbation equation for models with a nonzero cosmological constant are 
given in Edwards & Heath (1976), P80, Martel (1991), Lahav et al (1991), and Bildhauer, Buchcrt, & Kasai 
(1992) (see also Carroll, Press, & Turner 1992) For flat models, the solution is 

5/ l\ 1/2 } w^ 2 dw 
o 

and for critical models, 

5+ = \—^ ) J (l + u,)3/2(l- 2M ,)3 - (108) 



Since a is not an elementary function of t for these models, we cannot express 6+ as a function of t in closed 
form. Notice, however, that supercomoving variables still played an important role. If we had not been using 
the definitions of ao given by equations (85) and (90), the solutions (107) and (108) would show explicit 
dependences upon fi an d Qxo, as they normally do in comoving variables. 



8.2.4- Universe with a Radiation Background (n = 4) 

Solutions for models with a radiation background were derived by Meszaros (1974) and Groth & Peebles 
(1975). In supercomoving variables, the growing solution for flat models (that is, n = in eq. [83]) is 

6e" 2 * 



8.2.5. Universe with a Nonclumping, Nonrelativistic Component (n = 3) 

Models with a nonclumping, nonrelativistic component (n = 3) were discussed by Bond, Efstathiou, 
& Silk (1980) and Wasserman (1981). Such models are of interest in the context of models like the Cold+Hot 
Dark Matter model (CHDM) in which one component of matter, the Hot Dark Matter, has a much larger 
effective Jeans length than the others. In that case, while fluctuations are able to grow in the other com- 
ponents on scales smaller than this effective Jeans length in the HDM, they do not grow in the HDM. The 
n = 3 model, therefore, becomes a useful approximation for small wavelengths in the CHDM model. For flat 
models, the solution is 

6+ = a m = (-t)- m / 2 , (110) 



where 

1 

m= 4 



l + (l + 24f> ) 1/2 l . (Ill) 
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For non-flat models, the solutions arc not power laws. In the limit a/a <C 1, however, these solutions 
approach a power law, but with a different exponent, 



1 

m = — 

4 



1 + 



24O 



1/2' 



(112) 



Notice that the solutions (110)— (112) do not approach 5+ oc a at early time as they do in all models 
with n < 3. Hence, there is no preferred normalization for these solutions. 

Unlike in the previous cases, the solutions (110)— (112) explicitly depend upon the density parameter fi - 
In all the previous cases, the solutions were not power laws. Hence, there was a natural time scale in each 
solution, enabling us to rescale the time variable in order to make these solutions identical. But in the case 
of equation (110), the solution is a power law, and therefore has no natural time scale that could be used 
for rescaling. Hence, the solutions (110) for different values of fio are genuinely different, and cannot be 
transformed into one another by a rescaling of the variables. 



8.2.6. Universe with Domain Walls (n = 1) 

For models with domain walls (n = 1), we could not find any published solution in the literature. We 
shall therefore derive our own solutions of equation (105) for n = 1. As in §7.3, we consider two specific 
cases, flat models and critical models. 

First, consider as a trial solution the decaying solution for the Einstein-dc Sitter model (Oo = 1; 
Qxo = 0), that is, 5- = x~ z / 2 . We substitute this solution into equation (105), and get 

|(i -n - n X0 )x 1 / 2 + ^n xo (n - i)x b ' 2 - n = o . (113) 

This trial solution is the correct solution if both terms in equation (113) vanish. The first term vanishes for 
fiat models. The second term vanishes either for n X a = or n = 1. The former case corresponds to the 
Einstcin-de Sitter model, the latter to a model with domain walls. Hence the decaying mode is the same in 
a flat universe with domain walls as in a flat universe with no X component. Of course, this is true only if 
we express 6- as a function of a (or x = a/a ). The relationship between a and t is model-dependent, and 
so is the relationship between 5 and t. 

We thus have found the decaying mode for flat models with domain walls. To find the growing mode, 
we set 5 = 5- J f(x)dx and substitute in equation (105), with 1 — £! — Qxo = and n = 1. We get a first 
order equation for f(x), which is easily solved. The final solution is 

a 

r 5 -3/2 f W 3/2 dw 

s+ = r f Jo^w (114) 



For critical models with domain walls, we could not find any elementary solution for equation (106). 
However, we can express the growing solution as an infinite series. Since the solution must take the form 5+ cx 
a in the limit n<l (this is true for all models with n < 3), we define 

oo 
1=0 
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We substitute this solution into the perturbation equation, and solve for the coefficients. The solution is 

, 8 41(? + 1)q_ 1 + (-2; 2 + Z + 1)q_ 2 
co = l, * = -, Q = w — ) , l>2. (116) 

As in all the previous cases with n < 3, we normalize these solutions by imposing 6 + = t~ 2 at early time. 



9. THE COLLISIONLESS PANCAKE PROBLEM 

For a ID, planar density perturbation, it is possible to solve the fully nonlinear problem exactly, up 
to the moment of caustic formation, in the limit in which pressure can be neglected (Sunyaev & Zel'dovich 
1972; Doroshkievich, Ryabenkii, & Shandarin 1973). This is the problem which is sometimes described as 
the cosmological "pancake" problem (Zel'dovich 1970). Writing this solution in supercomoving variables 
serves to demonstrate the universal character of the solution, independent of cosmological model, as follows. 
We consider, for simplicity, a ID, primordial plane-wave density fluctuation as the initial condition, which 
is a solution of the linear perturbation equations discussed in the previous section. 

In the collisionless limit (p — 0), the fluid conservation equation in supercomoving form, for plane- 
symmetric systems, reduce to 

dv „di> dd> , ^ 

SU^fl-iY (ii 9) 



dx 2 \p 

We now set up initial conditions at time U by displacing each mass element from its unperturbed location 
q by a small amount as follows, 

ii = q- ii . (120) 
For the particular problem of the collisionless pancake, we choose initial conditions that are periodic, 

^ = -\sm2nkq, (121) 
2irk 

where k is the wavenumber of the perturbation. The corresponding initial density profile is given by 

~ Pl = 1 =r- ~^(l + 6 jC os27rfcq), (122) 

1 — bi cos 2irkq 

where the last equality is valid in the limit bi <C 1. The initial density contrast is therefore given by 
Si = bi cos 2nkq. Since the density contrast at early times grows according to the solution of the linear 
perturbation equations, the solution at early times must reduce to 

S(t) = b,^t\ cos 2irkq . (123) 
M*») 
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Our goal is to find a solution of the fluid equations (117)— (119) which reduces to equation (122) at early 
times. The solution is given by 

i(x,t) = ^lsm2wkq, (124) 
2nk 

where the amplitude b(t) is equal to bi at the initial time. We differentiate equation (120), and get 

~ = _db/dt sin27r ~ k ^ (125) 
27rfc 

We now substitute equation (125) into the continuity equation (117), and solve for the density, according to 

P 



1 — bcos2-Kkq 



(126) 



We then substitute equations (125) and (126) into equations (118) and (119). We are left with two unknowns, 
<f> and b(t). We eliminate <f>, and get 

fb 
dt- 



-^=6f n ab. (127) 



This equation is identical to equation (104). The solution is therefore given by 

b(t) = AS+(t) + BS-(t) . (128) 

By imposing the boundary condition (123), we get 

A = h/S+ik), (129) 
B = . (130) 

In supercomoving variables, the solutions for 5+ do not depend upon the cosmological parameters 
within each family, except for the case n = 3. The solution (124)-(126) describes the growth of a sinusoidal 
perturbation that forms a "pancake," that is, a caustic surface of infinite density, when the amplitude b(t) 
reaches unity. The solution is independent of the cosmological parameters within each family. When b{t ) 
reaches unity, this solution breaks down. At this point, strong shock waves form on both sides of the pancake, 
and the collisionless approximation (p = 0) is not valid anymore. Notice that the fluid equations are exactly 
the same in supercomoving variables as in noncosmological variables. This implies that the jump conditions 
in supercomoving variables do not depend upon the cosmological parameters. Since the solutions before 
shock formation do not depend upon these parameters either, we conclude that the postshock solution will 
also be independent of the cosmological parameters within each family. This means, in particular, that 
the pancake collapse problem, in the absence of non-adiabatic physical processes like radiative cooling, has 
one solution only, when expressed in supercomoving variables. Hence, for pancakes occurring in different 
cosmological models, the values and spatial variations of any fluid variable at any time are exactly the same 
as long as length is expressed in units of the pancake wavelength while time is expressed in terms of the 
cosmic scale factor a/a c , where a c is the scale factor at which b(i) = 1 and caustics are predicted to occur as 
described above, independent of the cosmological model. Approximate analytical solutions of the pancake 
collapse problem both before and after shock formation will be presented in a forthcoming paper (Shapiro, 
Struck-Marcell, & Martel 1998). 



10. THE ZEL'DOVICH APPROXIMATION 



Zel'dovich (1970) presented an analytical approximation for the growth of density perturbations, which 
extrapolates the linear solution into the nonlinear regime. This approximation is exact in the linear regime 
and, for locally ID motion, remains valid into the nonlinear regime, up to the formation of the first caustics. 
In this section, we derive the supercomoving form of the Zel'dovich approximation. 

Let r(t) be the proper position of a mass element at proper time t, and let r 4 = r(^) be the position 
of that mass element at some initial time U, in the absence of perturbation. It can be shown that, at early 
time, the proper position r is given by 



r(t) =a{t) q + <y+(t)p(q) 



(131) 



where q = r i/a(t j) , and p is a vector field that is function of q only (not t) , and characterizes the perturbation. 
The quantity dpj/dqu is sometimes referred to as the deformation tensor Djk- The time-independence of 
p implies that at early times mass elements are moving on straight lines. The Zel'dovich "approximation" 
consists of assuming that equation (131) is valid, not only at early times, but up to the point where mass 
element trajectories cross, at which point the density diverges. Writing equation (131) in supercomoving 
form yields 

f(t) = q + J+(t)p(q). (132) 

Notice that q and p were computed using equation (15) without the factor of a, since q and p in equation (131) 
are already expressed in comoving variables. In supercomoving variables, q = [a/a(ii)]fj, and the deformation 
tensor is unchanged, Djk — Djk- As in the case of the collisionless pancake problem, we get a solution whose 
dependence upon the cosmological model is entirely contained in the growing mode 5 + . Therefore, the 
solutions are identical within any given family. 



11. SUMMARY AND CONCLUSION 

In this paper, we described a set of dimensionless variables, which we call supercomoving variables, for 
the description of the state and evolution of matter in a cosmologically expanding Fricdmann universe in 
the Newtonian approximation (i.e. nonrelativistic gas on subhorizon scales). These supercomoving variables 
have the following properties: (1) Thermodynamic variables which depend in general upon space and time 
are stationary in the absence of perturbation, so that the effects of universal adiabatic cosmic expansion are 
accounted for implicitly and do not appear explicitly as sources of time-dependence, as long as the ratio 
of specific heats is 7 = 5/3. (2) Velocities reflect departures from uniform Hubble flow only and in the 
absence of peculiar gravity or pressure gradients, are stationary. (3) One advantage of these variables over 
standard comoving variables is that, for ideal gas with ratio of specific heats, 7 = 5/3, the fluid conservation 
equations in supercomoving variables are identical to the Newtonian fluid conservation equations for a 
noncosmological gas, while in standard comoving variables, extra terms appear to describe the effects of 
cosmology. (4) Similarly, the collisionless Boltzmann equation in supercomoving variables is identical to that 
for a nonrelativistic, noncosmological, collisionless gas. (5) The choice of nondimensionalizing units removes 
explicit dependence of the variables and the equations which describe their evolution including the Poisson 
equation for peculiar gravitational field, on the Hubble constant Hq. 

Our main objective in deriving supercomoving variables was to obtain a set of cosmological fluid equa- 
tions that resemble as closely as possible the fluid equations in noncosmological variables, and do not show 
any explicit dependences upon cosmological parameters such as H , Q , &xo, OT n. This was first achieved 
for matter-dominated models by Shandarin (S80) for adiabatic gas, with terms for nonadiabatic physical 
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processes added by Shapiro et al. (1983) and Shapiro & Struck-Marcell (1985). Here we have extended this 
approach to include models with a generalized, smooth background energy density (such as a cosmological 
constant) as well as additional physical processes, like viscosity and magnetic fields. The continuity equation, 
the momentum equation, and the equation of state have exactly the same form in noncosmological and su- 
percomoving variables. The cosmological drag term, that is present in the momentum equation in comoving 
variables, is absent in supcrcomoving variables. The energy equation has the same form in noncosmological 
and supercomoving variables, except for an additional drag term (which is also present in comoving vari- 
ables). However, unlike the case of standard comoving variables, that drag term vanishes for the physically 
interesting case 7 = 5/3. Only the Poisson equation remains very different from its noncosmological form. 
In supercomoving variables, as in comoving variables, the source of the gravitational potential is not the 
density, but rather the density contrast. 

The supercomoving variables reveal the existence of similarities among various cosmological models. 
This was shown in S80 for matter-dominated models. In comoving variables, the fluid equations, and 
their solutions, depend explicitly upon the density parameter fio- In supercomoving variables, all matter- 
dominated models fall into three categories, called families, defined by Clo < 1> f^o = 1, an d Sl > 1- Within 
each family, the fluid equations and their solutions are independent of £! - Therefore, using supercomoving 
variables effectively reduces the number of matter-dominated models from infinity to 3. The existence of these 
families has direct practical significance for the solution of problems involving gas dynamics and gravitational 
clustering in cosmology, as follows. Any solution of the differential equations in supercomoving variables 
obtained for a given initial or boundary value problem for a particular set of values of the cosmological 
parameters H and Slo for which those initial and boundary values have also been expressed in supercomoving 
variables immediately yields the solution for any other values of H and fl , within the same family as well, as 
long as the supercomoving initial and boundary values are also the same. This "Family" membership, then 
implies a "similarity" solution for initial and boundary value problems, where the supercomoving variables 
are the similarity variables. 

In this paper, we generalized supercomoving variables to cosmological models with two components, the 
usual nonrelativistic component, and a uniform component, called the X component, whose energy density 
varies as a~ n . This generic two-component model includes as special cases several physically motivated 
cosmological models, such as models with a nonzero cosmological constant, domain walls, cosmic strings, 
massive neutrinos, or a radiation background. Each of these model is a two-parameter model, defined by 
the density parameters f^o and fixo of the nonrelativistic and X components, respectively. 

We extended the concept of family introduced in S80 for matter-dominated models. We showed the 
existence of similarities among models with different values of f2 and ilxo- For any particular value of n, 
the various families arc defined by the value of the parameter k appearing in equation (71). Hence, each 
combination (n, n) corresponds to a particular family, and, as we showed, the fluid equations and their 
solutions do not depend upon the cosmological parameters f2o and Qxq within a particular family. While 
Shandarin (S80) reduced a one-parameter model (the parameter being Slo) to three no-parameter models, 
we are reducing two-parameter models (one for each value of n, the parameters being f2 and Qxo) to 
one-parameters models (the parameter being k). Therefore, as in the case of matter-dominated models, 
solutions obtained in supcrcomoving variables for a particular model (that is, for particular values of Qo 
and fijco) can be applied to other models as long as these models are in the same family. In this paper, wc 
focused our attention on 14 particular families which are physically interesting: (1) the Einstein-de Sitter 
model, (2) open, matter-dominated models, (3) closed, matter-dominated models, (4) flat models with a 
nonclumping background of nonrelativistic matter, such as would describe the case of massive neutrinos in the 
small wavelength limit, (5) open models with a nonrelativistic matter background, such as massive neutrinos, 
(6) closed models with a nonrelativistic matter background, such as massive neutrinos, (7) marginally bound 
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models with cosmic strings, [| (8) unbound models with cosmic strings, (9) bound models with cosmic 
strings, (10) flat models with a radiation background, (11) flat models with a nonzero cosmological constant, 
(12) critical models (i.e. asymptotically static) with a nonzero cosmological constant, (13) flat models with 
domain walls, and (14) critical models (i.e. asymptotically static) with domain walls. Table 4 summarizes 
the values of the parameters that define these various families. 

To illustrate the use of supercomoving variables and the value of the "family" membership they reveal, 
we have, for each of these families, described the linear growth of density fluctuations in supercomoving 
variables, the supercomoving Jeans length, and the nonlinear growth of ID planar density fluctuations - the 
cosmological pancake problem - prior to caustic formation. The solutions of these problems are independent 
of the particular values of the cosmological parameters f2o and Clx o within a given family, with the exception 
of the models with a nonrelativistic matter background such as massive neutrinos, for which the solutions 
must remain explicitly parametrized by Qq and flxo- The model with massive neutrinos differs from the 
other models in that the energy density of the two components both vary as a -3 . As a result, there is 
no characteristic time scale when the universe goes from being matter-dominated to being X component- 
dominated, or vice- versa. This absence of characteristic time scale implies that solutions for different models 
are genuinely different, and cannot be turned into one another by a scaling transformation, as for the other 
models. The only exception is the evolution of the scale factor (eq. [76]), which is made independent of Qq 
and tlxo by our introduction of the parameter /„ in the definition of the fiducial time i* (eq. [25]). By this 
choice, however, the form of the Poisson equation is different for n = 3 models. 

Finally, we have shown how the Navier-Stokes equation and the viscous energy equation are expressed in 
supercomoving variables, and presented supercomoving forms for the inviscid equation for vorticity evolution 
and the MHD equations, including the Biermann battery term. The latter two make clear how the vorticity 
and cosmic magnetic fields generated by large-scale structure formation, as discussed recently by Kulsrud et 
al. (1997), are related to each other and how they scale with characteristic time and length scales for cosmic 
structure formation in most cosmological models of interest. 



This work was supported by NASA Astrophysics Theory Grant Program Grant NAG5-2785 and NSF 
Grant ASC9504046. 



2 There is an important distinction, often forgotten, between the terms "closed" and "bound." The former 
refers to the curvature of the universe, and is determined by the sign of Qq + Qxo ~ 1; the latter refers to 
whether the universe will expand forever or eventually recollapse. For matter dominated models (and also 
models with massive neutrinos), closed models are always bound and open models are always unbound. But 
this is not true in general for other models. 
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A. ADDITIONAL PHYSICAL PROCESSES 

A.l. Heating, Cooling, and Thermal Conduction 

Shapiro ct al. (1983), and Shapiro & Struck-Marcell (1985) showed how to transform the terms for the 
additional physical processes of radiative cooling, external heating, and thermal conduction into superco- 
moving variables. This requires only a modification of the energy equation. In noncosmological variables, 
this equation becomes 

ds „ r-A V-S 

— + u- Ve = -^V-u+ , Al 

at p p p 

where V is the volume heating rate, A is the volume cooling rate, and S is the conductive flux. We define 
the supercomoving heating rate F, cooling rate A, and conductive flux S as 

f (A2) 



a 7 A 

A =y, (A3) 
S = ^ , (A4) 



where 



r* = A, = £, (A5) 
S* = p*vl. (A6) 

Equation (Al) becomes, in supercomoving variables, 

-= +v- Ve = -W(37- 5)e- ■ v+ — — . (A7) 

at p p p 

This equation has the same form as equation (Al), except for the drag term, which vanishes for 7 = 5/3. 



A. 2. Viscosity 

If the fluid has a finite viscosity, the momentum equation is replaced by the Navier-Stokes equation, 

£ + (u • V)u = - V* - ^ + ^V 2 u + ^±^V(V • u) , (A8) 
at p p p 

(Landau & Lifshitz 1987) where r\ and ( are the coefficients of viscosity. The presence of viscosity dissipates 
kinetic energy into heat, hence the energy equation must be modified accordingly. Its final form is 

§ + u.V £ = -?V.u+I^, (A9) 
at p p J orj 

where the stress tensor a' is defined by 

, fdui duj 2 du k \ du k 

a - = \-bv J + ^--f^) + c% ^ • (A10) 
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Wc define the following supercomoving coefficients of viscosity and stress tensor: 

V = ^, (All) 

T]* 

C = (A12) 



< - iw J + wr^) + Hd +6H )> (A13) 

where 

V* = (* = ^ • (A14) 
With these definitions, equations (A8) and (A9) become, in supercomoving variables, 

2 + (v • V)v = -n 5 + ?V 2 v + ^±^V(V ■ v) , (A15) 
ot p p p 

$ + v • W = -W(37 - 5)e - ?V • v + + ^W 2 . (A16) 

A.3. Vorticity 

The vorticity is defined by 

w = Vxu. (A17) 
In the absence of viscosity, the vorticity evolves according to 

— + V x (w x u) = -Vp x V( -J . (A18) 
To convert these equations to supercomoving form, we define 

a = ^ , (Ai9) 



where 

1 

w* = — 



(A20) 



Equations (A17) and (A18) become, in supercomoving variables, 

£ = Vxv, (A21) 

and 

— = + V x (u; x v) = -Vp x V - . (A22) 

at \pj 
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A.4. Magnetic Fields 

In the presence of magnetic fields, the momentum equation must be modified to take into account the 
magnetic tension and pressure. In the MHD approximation (where the fluid is assumed to have infinite 
conductivity), the momentum equation becomes 

<9u , „^ Vp (VxB)xB , 4 „ oN 
— + u • V u = -V* - — + ^ - 1 , A23 

ot p Airp 

where B is the magnetic field. The time evolution of the magnetic field is described by the following equation, 

— = V x (u x B) + -^-S7p e x Vn e . (A24) 

CjXi ^ 

This is a standard MHD equation, except for the last term, which we have added to take into account 
the possibility of magnetic field generation in an ionized gas by a process known as the Bicrmann battery 
(Bicrmann 1950). In this term, p e and n e are the electron pressure and number density, respectively. If we 
assume that the ionized fraction \ is uniform and the electron, ion, and neutral temperatures are equal, this 
equation reduces to 

dB = V x (u x B) - -^-Vp x vf -), (A25) 



dt e(l + X ) \p 

(Kulsrud et al. 1997), where toh is the hydrogen mass. To convert these equations to supercomoving 
variables, we consider two different definitions for the supercomoving magnetic field, 

-D* 

B 2 = (A27) 



where 



= P l /2 v* . (A28) 



Using the definition (A26), the momentum equation becomes, in supercomoving variables, 

ot p Airp 

which has the same form as equation (A23), and equation (A24) becomes 

r B! = V x (v x Bx) - — _JL_VpxV - . (A30) 



dt 2 v L > p i/2 n e(l + X ) \P 

Comparing equations (A24) and (A30), we see that the conversion to supercomoving variables introduces 
an "antidrag" term TiS/2 (which also appears in standard comoving variables, though with a different 
coefficient), and an extra factor of a 1 / 2 in the Biermann term. Using the second definition (A27), instead, 
the momentum equation becomes 

g + (v.V)v = -^-V0+ a ^ X ^ )x ^ , (A31) 

Ot p 47T/9 
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and equation (A24) becomes 



^ = V x (v x B 2 ) - ' VP x vf V] , (A32) 



0* pi /2 r,e(l + x) VP 

which has a form very similar to equation (A24). 

As noted recently by Kulsrud et al. (1997), there is a close correspondence between the inviscid equation 
for vorticity evolution, equation (A18) and the MHD equation for B in the presence of the Biermann battery 
term, equation (A25), especially when B is expressed in the same units as vorticity, by replacing B by the 
cyclotron frequency u> cyc = eH/runC. In the case in which both vorticity and magnetic field are initially 
zero, this leads to the remarkably simple result that 



1 + X 

We can write equation (A32) in terms of a> cyc defined by 



(A33) 



^cyc* 

where 



£ cyc = ^^, (A34) 



^cyc* - — . — ^* • (A35) 



This yields 



dti cvc ~ _ N 1 ~ ~fl 



^ - V x (v x £j cyc ) - -r^-Vp x V ( ^ ) . (A36) 



3t v w 1 + x r \P 

A comparison of equations (A22) and (A36) then replaces equation (A33) by 

U cyc = —— • 

1 + X 



(A37) 
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Table 1. SUPERCOMOVING AGES OF THE UNIVERSE (PRESENT AND ASYMPTOTIC) FOR 

FLAT MODELS WITH n = AND n = l. a 



n 




n = 






n = 1 








to 






to 




0.01 


4.6261 


0.20324 


0.21490 


9.9499 


0.27290 


0.28350 


0.02 


3.6593 


0.19630 


0.21490 


7.0000 


0.26558 


0.28350 


0.05 


2.6684 


0.18015 


0.21490 


4.3589 


0.24728 


0.28350 


0.10 


2.0801 


0.15834 


0.21490 


3.0000 


0.22080 


0.28350 


0.20 


1.5874 


0.11014 


0.21490 


2.0000 


0.17131 


0.28350 


0.50 


1.0000 


0.00000 


0.21490 


1.0000 


0.00000 


0.28350 


0.80 


0.6299 


-0.21178 


0.21490 


0.5000 


-0.33823 


0.28350 


0.90 


0.4807 


-0.38043 


0.21490 


0.3333 


-0.63233 


0.28350 


0.99 


0.2162 


-1.07573 


0.21490 


0.1005 


-2.02892 


0.28350 



a Present scale factor ao = a{io), where to = t(to), to = present 
age of the universe in proper time, and = t(t = oo), where t is 
proper time. We note that t(t = 0) = — oo in all cases. 
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Table 2. SELECTED FAMILIES OF INTEREST 



j_ eiinii y 


1 b 


n. 




Einstein-de Sitter model 


NA a 


NA 


Q = 1, n xo = 


Open, matter-dominated models 


NA 


NA 


Q < 1, Vxo = 


Closed, matter-dominated models 


NA 


NA 


O > 1, Qxo = 


Flat models with massive neutrinos 


3 


NA 


Qo + Qxo = 1 


Open models with massive neutrinos 


3 


NA 


n Q + n xo <i 


Closed models with massive neutrinos 


3 


NA 


n + n xo >i 


Marginally bound models with cosmic strings 


2 


NA 


O = 1, any Q xo 


Unbound models with cosmic strings 


2 


NA 


O < 1, any Q xo 


Bound models with cosmic strings 


2 


NA 


Cto > 1, any Q xo 


Flat models with a radiation background 


4 





flo + ^xo = 1 


Flat models with a nonzero cosmological constant 








flo + Qxo = 1 


Critical models with a nonzero cosmological constant 





-3/4V3 


Equation (68) 


Flat models with domain walls 


1 





Qo + = 1 


Critical models with domain walls 


1 


-2 


Equation (68) 



a For entries labeled "NA," this means that the parameter is "not applicable. 
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Figure Captions 



Fig. 1.— Evolution of Cosmological Models in Supercomoving Variables, (a) (Top panel) Cosmic Scale 
Factor a(t ) versus time t for various models. Solid curves (from left to right): open, flat, and closed, matter- 
dominated models. The solid dot indicates the location of the present for the flat model. Dotted curve: 
flat model with radiation. This curve merges with the solid curve for the flat matter-dominated model at 
late time. Short-dashed curves: critical models with a nonzero cosmological constant (lower curve) and with 
domain walls (upper curve). Long-dashed curves: flat models with a cosmological constant and with domain 
walls, (b) (Bottom panel) Supercomoving Hubble parameter H = a~ 1 da/dt versus t, for the same models 
as plotted in (a) 

Fig. 2. — Families of Cosmological Models (a) (Top left panel) Family membership curves for n = 0. The 
curves are isocontours of the family membership parameter n, as defined by equation (71), for particular 
values of k, as labeled. The dashed curves represent particular cases mentioned in the text. For values in 
the shaede ragion, the universe does not begin with a Big Bang, (b) (Top right panel) Enlargement of lower 
right corner of (a). Models located below the dashed curve are bound, others are unbound, (c) (Middle left 
panel) Same as (a), except for models with n = 1. (d) (Middle right panel) Family Memberships for n = 2 
models. There are 3 families, represented by the dashed line, and the 2 half-planes, (e) (Bottom left panel) 
Same as (d) except for models with n = 3. (f) (Bottom right) Same as (a) except for models with n = 4. 
Shaded area in (a) and (c) represents combination of parameters incompatible with the existence of a Big 
Bang. 
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